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ABSTRACT 



The beamed radio emission from relativistic plasma (particles or bunches), 
constrained to move along the curved trajectories, occurs in the direction of ve- 
locity. We have generalized the coherent curvature radiation model to include the 
detailed geometry of the emission region in pulsar magnetosphere, and deduced 
the polarization state in terms of Stokes parameters. By considering both the 
uniform and modulated emissions, we have simulated a few typical pulse pro- 
files. The antisymmetric type of circular polarization survives only when there is 
modulation or discrete distribution in the emitting sources. Our model predicts 
a correlation between the polarization angle swing and sign reversal of circular 
polarization as a geometric property of the emission process. 

Subject headings: polarization - pulsars: general — radiation mechanisms: 
non-thermal 
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1. Introduction 

Pulsars are highly magnetized with predominantly dipolar field structure. The rotating 
magnetic field produces a strong induced electric field that accelerates charged particles off 
the surface of the star into a magnetosphere consisting of predominantly dipolar magnetic 
field and corotating relativistic pair plasma. Pulsar radio emission models assume that 
radiation emitted tangentially to the field lines on which plasma is moving. The polarization 
state of the emitted radiation is more or less determined by the structure of magnetic 
field at the emission spot. In the general framework of models in which the radio power 
is curvature radiation emitted by charge bunches constrained to follow field lines, the 
linear polarization is intrinsic to the emission mechanism and is, furthermore, a purely 
geometric property. Several pulsar researchers have shown that the properties such as the 
polarization angle swing can be explained within the framework of curvature radiation 
(e.g., Radhakrishnan & Cooke 1969; Sturrock 1971; Ruderman & Sutherland 1975; Lyne & 
Manchester 1988; Rankin 1990, 1993; Blaskiewicz et al. 1991). 

The radio emission from particle bunches is highly polarized, and the radiation received 
by distant observer will be less polarized due to the incoherent superposition of emissions 
from different magnetic field lines (Gil & Rudnicki 1985). Gil (1986) has argued for the 
connection between pulsar emission beams and polarization modes and suggested that 
out of two orthogonal polarization modes one corresponds to core emission and other to 
the conal emissions. They are highly linearly polarized and the observed depolarization 
is due to superposition of modes at any instant (Gil 1987). By considering a charged 
particle moving along the curved trajectory (circular) confined to the xz-plane, Gil and 
Snakowski (1990a) have deduced the polarization state of the emitted radiation, and shown 
the creation of antisymmetric circular polarization in curvature radiation. By introducing a 
phase, as a propagation effect, the difference between the components of radiation electric 
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field in the directions parallel and perpendicular to the plane of particle trajectory, Gil and 
Snakowski (1990b) have developed a model to explain the depolarization and polarization 
angle deviations in subpulses and micropulses. Gil, Kijak and Zycki (1993) have modeled 
the single pulse polarization characteristics of pulsar radiation, and demonstrated that 
the deviations of the single pulse position angle from the average are caused by both 
propagation and geometrical effects. Mitra, Gil and Melikdze (2009) by analysing the 
strong single pulses with highly polarized subpulses from a set of pulsars, have given a very 
conclusive arguments in favor of the coherent curvature radiation mechanism as the pulsar 
radio emission mechanism. 

By analyzing the average pulse profiles, Radhakrishnan and Rankin (1990) have 
identified two most probable types of circular polarizations, namely, antisymmetric, where 
the circular polarization changes sense near the core region, and symmetric, where the 
circular polarization remains with same sense. They found that antisymmetric circular 
polarization is correlated with the polarization angle swing, and speculate it to be a 
geometric property of the emission mechanism. Han et al. (1998), by considering the 
published mean profiles, found a correlation between the sense of circular polarization and 
polarization angle swing in conal double profiles, and no significant correlation for core 
components. Further, You and Han (2006) have reconfirmed these investigations with a 
larger data. However, Cordes et al. (1978) were the first to point out an association between 
the position angle of the linear polarization and the handedness of the circular polarization. 

There are two types of claims for the origin of circular polarization: intrinsic to 
the emission mechanism (e.g., Michel 1987; Gil & Snakowski 1990a, b; Radhakrishnan 
& Rankin 1990; Gangadhara 1997) or generated by the propagation effects (e.g., Cheng 

& Ruderman 1979). Cheng and Ruderman (1979) have suggested that the expected 
asymmetry between the positively and negatively charged components of the magnetoactive 
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plasma in the far magnetosphere of pulsars will convert linear polarization to circular 
polarization. Radhakrishnan and Rankin (1990) have suggested that the propagation 
origin of antisymmetric circular polarization is very unlikely but the symmetric circular 
polarization appears to be possible. On the other hand, Kazbegi, Machabeh and Mehkidze 
(1991, 1992) have argued that the cyclotron instability, rather than the propagation effect, 
is responsible for the circular polarization of pulsars. Lyubarskii and Petrova (1999) 
considered that the rotation of the magnetosphere gives rise to wave mode coupling in the 
polarization-limiting region, which can result in circular polarization in linearly polarized 
normal waves. Melrose and Luo (2004) discussed possible circular polarization induced by 
intrinsically relativistic effects of pulsar plasma. Melrose (2003) reviewed the properties 
of intrinsic circular polarization and circular polarization due to cyclotron absorption and 
presented a plausible explanation of circular polarization in terms of propagation effects in 
an inhomogeneous birefringent plasma. In the multifrequency simultaneous observations 
we do find the variations in the single pulse polarization, which may be attributed to the 
propagation effects (Karastergiou et al. 2001, 2002; Karastergiou, Johnston & Kramer 
2003). 

The correlation between the antisymmetric circular polarization and the polarization 
angle swing is a geometric property of the emission processes (Radhakrishnan and Rankin 
1990). By carefully modeling the polarization state of the radiation in terms of Stokes 
parameters, it is possible to construct the geometry of emission region at multifrequencies. 
So far, in the purview of curvature radiation only the polarization angle has been modeled 
(Radhakrishnan & Cooke 1969; Komesaroff 1970) and attempted to fit it with the average 
radio profile data (e.g., Lyne & Manchester 1988). Instead of circular trajectories, it 
is very important to consider the actual dipolar magnetic field lines, whose curvature 
radii vary as a function of altitude, as the radio emission in pulsars is expected to come 
from the range of altitude (e.g., Gangadhara & Gupta 2001; Gupta & Gangadhara 2003; 
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Krzeszowski et al. 2009). In this paper, we develop a three-dimensional (3D) model for 
curvature radiation by relativistic sources accelerated along the dipolar magnetic field lines. 
We consider the actual dipolar magnetic field lines (not the circles) in a slowly rotating 
(non-rotating) magnetosphere such that the rotation effects can be ignored. The relativistic 
plasma (bunch, i.e., a point- like huge charge) moving along the dipolar magnetic field lines 
emits curvature radiation. We show that our model reproduces polarization angle swing 
of Radhakrishnan and Cooke (1969), and predicts that the correlation of antisymmetric 
circular polarization and polarization angle swing is a geometric property of the emission 
process. Our model is aimed at re-examining the intrinsic polarization properties of the 
vacuum single-particle cuvature radiation, and planed to consider the propagation effects 
separately in the subsequent works. We derive electric fields of the radiation field in 
section 2 and construct the Stokes parameters of the radiation field in section 3. A few 
model (simulated) profiles are presented in section 4 depicting the correlation between the 
antisymmetric circular polarization and polarization angle swing in the different cases of 
viewing geometry parameters. 

2. Electric field of curvature radiation 

Consider a magnetosphere having dipole magnetic field with an axis m inclined 
by an angle a with respect to the rotation axis VL (see Figure [1]). We assume that the 
magnetosphere is stationary or slowly rotating such that the rotation effects are negligible. 
The relativistic pair plasma, generated by induced electric field followed by pair creation, 
constrained to move along the curved dipolar magnetic field lines. The high brightness 
temperature of pulsar indicates coherency of the pulsar radiation, which in tern forces one 
to postulate the existence of charged bunches. The formation of bunches in the form of 
solitons has been proposed (e.g., Cheng & Ruderman 1979; Melikidze & Patarya 1980, 1984) 
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and questioned (e.g., Melrose 1992). Gil, Lyubarsky and Melikidze (2004) have generalized 
the soliton model by including formation and propagation of the coherent radiation in the 
magnetospheric plasma along magnetic field lines. Their results strongly support coherent 
curvature radiation by the spark-associated solitons as a plausible mechanism of pulsar 
radio emission. Following these views, we assume that the plasma in the form of bunches 
moves along the open field lines of the pulsar magnetosphere. 

Consider the source S moving along the magnetic field line C, and experiencing 
acceleration (a) in the direction of curvature vector of the field line. We assume that the 
source to be a bunch, which is nothing more than a point-like huge charge. In Cartesian 
coordinates, the position vector of a bunch moving along the dipolar magnetic field line is 
given by (see Equation (2) in Gangadhara 2004, hereafter G04) 

r = Tg sin^ 0{cos cos 0' sin a + sin ^(cosa cos (/)cos(/)' — sin 0sin0'), 
cos (f)' sin 6 sin -|- sin 0' (cos 9 sin a + cos a cos (j) sin 9) , 
cos a cos 9 — cos sin a sin 9} , (1) 

where r^. is the field line constant, and the angles 9 and are the magnetic colatitude and 
azimuth, respectively. Next, 0' is the rotation phase and a is the inclination angle of the 
magnetic axis. Equation ([T]) describes the dipolar magnetic field lines presented in Figure 1. 
Then the velocity of bunch is given by 

dr {dr\ fd9\ fd9\ , 

where b = dr/d9 is the magnetic field line tangent. Consider the magnetic axis: 

m = {sina COS0', sin a sin 0', cos a}. (3) 

Due to curvature in the field lines, the plasma bunch, a point-like huge charge, 
collectively radiates relativistically beamed radiation in the direction of velocity v. The 
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velocity v is parallel to the tangent b of the field line. To receive the beamed emission, the 
observer's line of sight (n) must align with v within the beaming angle I/7, where 7 is 
the Lorentz factor of the bunch. In other words a distant observer at P receives beamed 
emission only when n ■ v = cosr ~ 1 for r ^ I/7, where v = v/\v\. Let s be the arc 
length of the field line then ds = \b\d9, where \b\ = (rg/ a/2) sin 6^/5 + 3 cos(20) , and the 
magnitude of velocity v = ds/dt = kc, where the parameter k specifies the speed of bunch 
as a fraction of the speed of light c. Hence, we have 

V = Kch , (4) 

where 

h = b/\b\ = {cosrcos0'sina + sinr(cosacos0cos0' — sin0sin0'), 
cos (f)' sin T sin (j) + sin 0' (cos r sin a + cos a cos (f) sin r) , 
cos a cos r — cos sin a sin r} , (5) 

and T is the angle between m and b. In terms of polar angle 6, the angle r given by 

sinr 3sin(2^) 

tanr= = . ^ ^ ^ ' , (6) 

cosr l + 3cos(26'j 

where 

1 + 3cos(2e) 



cos T = h ■ rh 



sin r = {m X b) ■ 



v/lO + 6cos(20)' 
3sin(20) 



^10 + 6cos(2^) 
and 

e<^ = {— cosa sin0 COS0' — cos sin0', cos0cos0' — cosQ;sin0sin0', sin a sin 0} (7) 
is the bi-normal to the field line. We solve equation ([H]) for 6, and obtain 

cos(26') = - (cos r VS + cos^ r — sin^ r) . (8) 
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Hence, from equation it is clear that to receive the radiation emitted in the direction of 
tangent b the hne sight hne n must hne up with it. So, by solving n ■ b = 1 oi n x b = 0, we 
can identify the tangent b, which aligns with n, and hence find the field line curvature and 
the coordinates (6', 0) of the emission spot (see Equations (4), (9), and (11) in G04). Next, 
the acceleration of the bunch is given by 

dv (kcY db 



a 



(kc) 



dt \b\ de 

where k = {l/\b\)db/d6 is the curvature (normal) of the field line. Then the radius of 
curvature of field line is given by 



(9) 



P 



1 

W\ 



2 - 



Therefore, using k = k/p, we can write 



3{3 + cos(2^)} 



[KC 



\b\ . 



(10) 



a 



P 



(11) 



where 



k = { (cos a cos (p cos (j)' — sin cj) sin 0') cos r — cos cp' sin a sin r, 
(cos 0' sin + cos a cos sin 0') cos r — sin a sin 0' sin r, 
— cos sin a cos r — cos a sin r} . 
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The relativistic bunch, i.e., point-like huge charge q that collectively emits curvature 
radiation as it accelerates along the curved trajectory C (see Figure 1). Then the electric 
field of the radiation at the observation point P is given by (Jackson 1975): 

n X [{h — f3) X /3] 



E{r,t) 



Re 



(13) 



ret 



where ^ = 1 — /3 • i? is the distance from the radiating region to the observer, f3 = v/c is 
the velocity, and = a/c is the acceleration of the bunch. 
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The radiation emitted by a relativistic bunch has a broad spectrum, and it can be 

estimated by taking Fourier transformation of the electric field of radiation: 

+00 

E{r,u) = ^ ! E{r,t)e"^'dt. (14) 



/27r 

In equation f|T3|) . ret means evaluated at the retarded time t' + R{t')/c = t. By changing the 
variable of integration from t to t', we obtain 



~oo 



E 



+00 . 



where we have used dt = ^ dt' . When the observation point is far away from the region 
of space where the acceleration occurs, the propagation vector or the sight line n can be 
taken to be constant in time. Furthermore, the distance R{t') can be approximated as 
R{t') ^ Rq — h ■ r{t'), where Rq is the distance between the origin O and the observation 
point P, and r{t') is the position of the bunch relative to O. 

Since bunches move with velocity k c along the dipolar field lines, over the incremental 
time dt the distance (arc length) covered is ds = ncdt = \b\d6. Therefore, we have 

t = 1- j \b\d9 = y"sin0y5 + 3cos(2^) dO . (16) 

By choosing t = at 6* = 0, we obtain 

12 + V31og (lA + sVfj - 3v^lO + 6cos(20) cos(^) - 

2^3 log (y6cos{9) + v/5 + 3cos(2e)) . (17) 

By assuming k ~ 1, in Figure [2], we plotted t as a function of 6 for different r^. It shows 
time t increases much faster at larger than at lower. This is due to the fact that for a 
given range of 6 the arc length of the field line becomes larger at higher r^. 

Then equation ( 1T5|) becomes 

+00 

Eir,.) . [ |,|^x[(n-/3)x/3]^... 



t 



12kc 



-00 
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where the expression for t is given by equation (ITTI) . Note that the prime on the time 
variable t has been omitted for brevity. The integration limits have been extended to ±00 
for mathematical convenience, as the integrand vanishes for — ^o| > 1/7- rotation 
phase 0', there exists a magnetic colatitude 60 and a magnetic azimuth 0o at which the field 
line tangent b exactly align with n, i.e., bo ■ h = 1 and r = F, where F is the half-opening 
angle of the pulsar emission beam centered on m. The expressions for ^0 and 0o are given 
in G04 (see Equations (9) and (11)). 

The polarization state of the emitted radiation can be determined using E{uj) with the 
known r{t), (3 and /3. Since the integral in equation f[T5]) has to be computed over the path 
of particle the line of sight h can be chosen without loss of generality, to lie in the a:2;-plane: 

n = (sin^, 0, cosC) , (19) 

where ( = a + a is the angle between h and Cl, and a is the closest impact angle of h with 
respect to rh. 

Let 



A = -|b| '^"^^'7^^^"^J . (20) 



1 n X [{n-f3) X f3] 

KC ^ 

By substituting for acceleration /3 = a/c from equation ([9]), we can reduce it to 

A = {A., Ay, = ^^f^ ^ , (21) 

where N = Kdb/dd= djS/dd. Using the expression (3 = v/c from equation (jlj) and series 
expanding A in power of 9 about 9q we obtain 

= A^o + A^,{e-eo) + AAe-eof + AAe-eof + o[{e-eof], 
Ay = A^o + Ay,{e-eo) + A^^{e-eof + A^^{e-eof + o[{e-eof], 
A = A,o + A,^{e-eo) + A,2{e-eof + A,:,{e-eof + o[{e-eoY], (22) 

where Axi, Ayi and A^i with i = 0, 1, 2, 3 are the series expansion coefficients, and their 
expressions are given in Appendix A. 



The scalar product between h and r is given by 



n 



■ r = Tf, sin^ 



cos a (cos 9 cos C, + cos cos 0' sin 9 sin (^) — cos C, cos sin a sin ^ + 



sin C (cos 6* cos 0' sin a — sin 6' sin sin 0') . (23) 

Next, substituting the expressions of t and n ■ r into the argument of exponential in 
equation (|T8l) . and series expanding in powers of 9 about 9o we obtain 

(71 ' T* \ 
t-^j =co + ci{9~9o) + C2{9 - 9of + c^{9 - 9of + 0[{9 - 9^f] , (24) 

where Cq, Ci, C2, and C3 are the series expansion coefficients, and their expressions are given 
in Appendix A. 

Now, by substituting the expressions of equations fl22|) and into equation ffTSl) . we 

obtain the components of E{uj) = {£'^(0;), Ey{u), Ez{uj)} : 

+00 

—00 
+00 

^y(w) = (AyO + /i + Ay2/i' + 3 /i^ ) 6^ ^^^^^ ^ ^ /i , 

^z(c^) = ^oy"(^zo + ^zi/i + ^z2/i' + ^3/i')e^^'=^'^+^^'^'+'=^^')rf/i, (25) 

—00 

where n = 9 — 9q and 

Now by substituting the integral solutions 5*0, Si, S2 and S^, given in Appendix B, into 
equation (l25ll we obtain 

Ex(w) = £^o(^xO^O + ^xl^l + ^x2^2 + ^xS^s) , 
Ey(w) = Eo(v4yo5o + ^yl5l + Ay252 + ^yS^s) , 

E,{uj) = Eo{A,oSo + A,iSi + A,2S2 + A,^S^). (26) 
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To find the polarization angle of radiation field E, we need to specify two reference 
directions perpendicular to the sight line h. One could be the projected spin axis on the 
plane of the sky: = (— cosC, 0, sin^), and then the other direction is specified by 
e_L = ey X n = y, where y is a unit vector parallel to the y-axis. Then the components of E 
in the directions ey and e± are given by 

E\\ = e\\ ■ E = — cos( + sin ( E^ , 

E± = e±-E = Ey . (27) 

At any rotation phase (f)', the observer receives the radiation from all those field lines 
whose tangents lie within the angle I/7 with respect to the sight line n. Let rj be the angle 
between the b and n, then cost] = b ■ n, and the maximum value of is I/7. Therefore, at 
= 00 we solve cos(l/ •y) = b - h for r, and find the allowed range (F — I/7) < r < (F + I/7) 
ofror— l/7<?7<l/7ofr7, which in turn allows to one find the range of 6 with the help 
of equation ([8]). Next, for any given rj within its range, we find by solving cos rj = b ■ n. It 
gives (00 — 50) < < (00 + (50), where 



, . . _ sinF[cos(l/7) csc(F + 77) - cosFcot(F + r])] 
cos(o0j — - — -— — — . 2 /- ■ 2 I, ■ 

(cos(, sma — cosa COS0 sm(^)^ + sm Qsm ^ 



Hence by knowing the ranges of 6 and at any given 0', we can estimate the contributions 
to E from all those field lines, whose tangents lie within the angle I/7 with respect to 
n. In Figure [3l we have plotted those regions at three phases; 0' = —30°, 0° and 30° 
using a = 10°, (3 = 5°, and 7 = 400. Note that at the center of each region b exactly 
aligns with the sight line, i.e., b ■ fi = 1. Further, in Figure HI we have plotted them for 
— 180° < 0' < 180° with a step of 5° between the successive regions. We observe that the 
range of 6 stays nearly constant (or decreasing negligibly) whereas that of gets narrower 
with respect to the increasing |0'|. 
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3. Polarization of radiation field 

To understand the pulsar radio emission, we must model all the Stokes parameters (/, 
Q, U and V — a set of parameters used to specify the phase and polarization of radiation), 
and compare with observations. They have been found to offer a very convenient method 
for establishing the association between the polarization state of observed radiation and the 
geometry of the emitting region. They are defined as follows: 

/ = + E^El , Q = EiiE^ - E±El , U = 2 RelE^Ej] , V = 2 lm[ElE±] . (29) 

The parameter / defines the total intensity, Q and U jointly define the linear polarization 
and it's position angle, and V describes the circular polarization. 



3.1. Addition of Stokes parameters 

Let Wj be the energy radiated coherently per unit solid angle per unit frequency 
interval per particle bunch (Jackson 1975), then 

= (30) 

Since the Stokes parameter / = E\\E\\* + E±E±* = E ■ E* = we can rewrite 

equation fl3Ql) as 

^ (31) 
Similarly, we can express Q, U and V as 

27r d'^Wq 



Q 



c Rq duj dVL ' 



^ 27r d'^Wt 



c -Rq du dQ ' 



^ 2r^fWy 

cRldudn ■ ^ ' 
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The net emission, which the observer receives along n, will have contributions from the 
neighboring field lines, whose tangents are within the angle I/7 with respect to n. Hence 
the radiation received at any given phase is the net contribution from a small tube of field 
lines having an angular width of about 2/7. Thus the radiation in the direction of n should 
be integrated over a sohd angle dQ — sin 9 dO dcf). We choose hmits on the angles and 9 
such that the integration over them will cover the solid angular region (beaming region) of 
radial width I/7 around h. Since and (j) are orthogonal, choosing them as the variables 
of integration is justified. We assume (i) the width of bunch r^o is much smaller than the 
wavelength A of the radio waves, so that the radiation emitted by a bunch is coherent, and 
(ii) the bunches, within the beaming region, are closely spaced, so that the net emission 
becomes smooth and continuous. 

Consider a bunch having 7 ~ 400 emitting radio waves at frequency u = 600 MHz at 
an altitude of about 400 km. Note that these values are closer to those estimated in G04 in 
the case of PSR B0329+54. Then the angular width of the beaming region corresponding 
to 2/7 is ~ 0.3°, which corresponds to width of ~ 2 km at an altitude of 400 km. For 
coherence to be effective the bunch width wq < A. Therefore, we choose Wq < 50 cm for 
A ~ 50 cm. Since these values of wq are much smaller than the width of the beaming region 
(~ 2 Km), the Stokes parameters can be integrated as continuous functions of 9 and (p. 

Let Is be the resultant Stokes intensity parameter then 

7s = y Idn 

'No- 
where 9q and 0o are the magnetic colatitude and azimuth of the sight line h. Similarly, for 
other Stokes parameters, we have 

Qs = / Qsm9d9d(j), 



I I Ism9d9d(t), (33) 
Jeo-se J<t>o-6<t> 
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'00-50 J(t)o-5<p 

V, = / Vsinededcp. (34) 



Then the hnear polarization is given by 



= VOT+U!, (35) 



and the corresponding polarization angle is 

1 



V's = J tan-i ] . (36) 



4. Simulation of pulse profiles 

The emission in spin-powered pulsars is mostly of non-thermal origin. If the radiation 
field E from different sources does not bear any phase relation then they are expected to 
be incoherently superposed on the observation point. On the other hand, if there is a phase 
relation then they are coherently superposed. From the observational point of view both 
the cases are important. 

By considering the relativistic pair plasma with 7 = 400 accelerated along the dipolar 
field lines of a pulsar with period P = 1 s, we computed the polarization parameters and 
plotted them in Figure [51 It shows a stronger emission near the meridional plane, where the 
beaming region is broader (see Figure [3]) and the radius of curvature p goes to a minimum. 
The profile of linear polarization Lg resembles the intensity profile, except for its lower 
magnitude due to the incoherent addition. To describe the behaviors of circular polarization 
and polarization angle ips, we define the symbols: '— /+' for transition of the right hand 
circular to left hand circular, '+/— ' for left hand circular to right hand circular, 'cw' for 
clockwise rotation of the polarization angle, and 'ccw' for counter clockwise rotation. 
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Since the circular polarization changes sign as — /+ or +/— as the sight line cuts 
across the field line, the net circular polarization goes to zero in an uniform emission due to 
the addition with opposite signs. The polarization angle swings reproduced in Figure O are 
consistent with the rotating vector model of Radhakrishnan and Cooke (1969). In the case 
of positive sight line impact parameter {a = 5°), the polarization angle swing is ccw as the 
slope d4's/d(p' > while in the negative case (a = —5°), it is cw as dips/dcp' < 0. 



4.1. Modulation of radio emission 

Pulsar radio emission is believed to come from mostly open magnetic field lines, whose 
foot points define the polar cap. The shape of pulsar profiles indicates that the entire polar 
cap does not radiate, only some selected regions radiate, which may be organized into 
a central core emission and coaxial conal emissions, which has an overwhelming support 
from observations (e.g., Rankin 1990, 1993). Hence the radiating region above the polar 
is believed to have a central column of emission (core) and a few coaxial conal regions 
of emission (cones) (e.g., Gil & Krawczyk 1997; Gangadhara & Gupta 2001; Gupta & 
Gangadhara 2003; Dyks, Rudak & Harding 2004). 



4- 1.1. Modulating function 

It is well known that the components of a pulsar profile can be decomposed into 
individual Gaussians by fitting one with each of the subpulse component. For example 
the components in the pulse profile of PSR 1706-16 and PSR 2351+61 are fitted with 
appropriate Gaussians by Kramer et al. (1994). When the line-of-sight crosses the 
emission region, it encounters a pattern in intensity due to Gaussian modulation in the 
azimuthal direction. Because of the Gaussian modulation in the azimuthal direction, the 
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intensity becomes nonuniform in the polar directions too. These arguments indicate that 
a Gaussian-hke intensity modulation exists in the polar directions too. So, we assume that 
the emission region of a pulse component has an intensity modulation both in the azimuthal 
directions. Hence we define a modulation function / for a pulse component as 



fi^A) = fo exp 



2 



(37) 



where 0p is the peak location of the Gaussian function and /o is the amplitude. If is the 
full width at half-maximum (FWHM) then = w^/{2y/h[2). 

Taking into account of modulation, equations fl5^ -f l51|) can be written as 

Is = / flsmeded(f) 



Jeo-se J(f>Q~5(i> 
Qs = / fQsmeded(j), 

J 00-50 Jipo-S^ 

f/s = / / f U sine ded(f) , 

J 00-50 Jipo-S^ 
1-00+50 r^o+5^ 

Vs = / fV sine ded(j) . (38) 

J 00-50 J(t>o-5<l) 

Using a Gaussian with peak located at the meridional plane (0' = 0°), we have 
computed the pulse profiles in the two cases of impact parameter (a) and inclination angles 
(a) and plotted in Figures E] and [3 We observe that the profile in the case of negative 
a is broader than the positive case. This difference is due to the projection of emission 
region on to the equatorial plane of the pulsar. In the case of positive a, the polarization 
angle Xs swing is ccw and the sign change of Vs is —/ + with respect to 0', while in the case 
of negative a the Xs swing is cw and the sign change of K is +/ — . Hence we find that 
the polarization angle swing is correlated with the circular polarization sign reversal. This 
correlation is invariant with respect the stellar spin directions. 

The mean pulsar profiles often found to consist of odd number of multi-components 
or subpulses. Many of the works on pulsar profiles (e.g., Rankin 1990, 1993; Mitra 
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& Deshpande 1999) have proposed that the pulsar emission beam has a nested conal 
structure. To investigate the polarization of sub-pulses in such profiles, we have reproduced 
a five component profile by considering three Gaussians in Figure [HI and five Gaussians 
in Figures |9] and [TOl The central component is presumed to be a core, and the other 
components are symmetrically located on either side of the core forming the cones. We find 
across each component that circular polarization changes the sign and is correlated with the 
polarization angle swing. We also observe that the circular polarization of the outermost 
components is weaker compared to that of inner ones, which is quite clear in the case of 
large inclination angles. This is due to the fact that the sight line crosses the field lines in 
the almost edge-on position in the case of outermost components. The small distortions in 
the polarization angle curve are due to modulation. 

5. Discussion 

Observed pulsar radio luminosities together with the small source size imply 
extraordinarily high brightness temperatures, i.e., as high as 10^^ K. The incoherent sum 
of a single-particle curvature radiation is not enough to explain the very high brightness 
temperature of pulsar radio emission, therefore, one is forced to postulate the existence 
of charged bunches. To avoid implausibly high particle densities and energies, coherent 
radiation processes are invoked. Pacini and Rees (1970), and Sturrock (1971) among others 
were quick to point out that the observed coherence may be due to bunching of particles in 
the emission region of the magnetosphere. The problem of bunch formation has been known 
for many decades, and it has already been addressed by many authors (e.g., Karpman et al. 
1975; Ruderman & Sutherland 1975; Cox 1979; Asseo, Pelletier & Sol 1990; Gil, Lyubarsky 
& Melikidze 2004). The natural mechanism for the formation of charged bunches was first 
proposed by Karpman et al. (1975). They have argued that the modulational instability in 
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the turbulent plasma generates charged solitons, provided that species of different charges 
have different masses. One should mention here that to explain coherent radio emission we 
do not necessarily need stable solitons but only large-scale (as compared with the Langmuir 
wavelength) charge density fluctuations. Gil, Lyubarsky and Mehkidze (2004) generahzed 
the soliton model by including formation and propagation of the coherent radiation in 
the magnetospheric plasma. However, it is not easy to form such charge bunches (see 
Melrose 1992 for a review). Further, Michel (1991) has pointed out that the pair-production 
discharge mechanism originally applied to pulsars by Sturrock automatically produces 
dense bunches that can produce coherence at radio frequencies with sufficient intensity to 
simulate pulsar action. If the bunches of plasma particles with sizes much smaller than a 
wavelength of radiation exist then the net radiation field £^(a;) fa NEo{uj), where N is the 
number of charges present in the bunch. Hence the total radiation field due to a bunch of 
particles is equal to the vector sum of the fields radiated by each charge. 

In the general framework of models including ours, in which the radio power is 
curvature radiation emitted by charge bunches constrained to follow the field lines, the 
linear polarization is intrinsic to the emission mechanism and is, furthermore, a purely 
geometric property. In this direction the recently acheieved observational results and the 
model predications based on them by Mitra, Gil and Mehkidze (2009) become very relevant. 
They find that the polarization angle of linear polarization in subpulses follow closely the 
mean polarization angle curve at the corresponding profile components and argue that their 
findings favor coherent curvature radiation over maser mechanism as the observed emission. 

In an actual case, it is the combination of both incoherent and coherent superpositions 
determining the polarization state of the observed emission. Though the emissions from a 
single bunch is highly polarized, the radiation received by a distant observer will be less 
polarized, as the radiation from many bunches is incoherently superposed (Gil & Rudnicki 
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1985). Also, the degree of polarization is found, depending on the time resolution chosen 
in the observation (Gangadhara et al. 1999). Circular polarization is generally strongest in 
the central regions of a profile, but is by no means confined to these regions. It has been 
detected from conal components of many pulsars, for example, conal-double pulsars and 
found to be highly correlated with the polarization angle swing (You & Han 2006). In most 
of the cases the sense reversal of circular polarization is nearly independent of frequency, 
suggesting that the circular polarization does not arise from propagation or plasma effects 
(Michel 1987; Radhakrishnan & Rankin 1990). Radhakrishnan and Rankin (1990) have 
argued that the circular polarization is intrinsically antisymmetric type and correlated with 
the polarization angle swing. The antisymmetric circular polarization of curvature radiation 
becomes significant if there are gradients in the emissivity over angular scales comparable 
with the emission cone of single charge. Their results are consistent with the predictions 
of our model (Figures [6] - [TOj) and strongly suggest that the correlation of antisymmetric 
circular polarization with the polarization angle swing is a geometric property of emission 
mechanism. Since our model deals with steady flow of relativistic plasma bunches along 
dipolar field lines, it is relevent only for average profiles, and reflects the results which are 
more of geometric dependent. We have not considered any fluctuations or instability in the 
plasma flow. Hence it can not reproduce the behaviors of single pulses. 

By adopting the antenna mechanism, Buschauer and Benford (1976) have derived a 
new formalism for the relativistic curvature radiation. However, the treatment given does 
not include the detailed geometry of dipolar field lines and the estimation of polarization, 
particularly circular, as we have considered in our model. Since radiation from many 
bunches is superposed on any given pulse longitude, the circular polarization of different 
signs and magnitudes is added. The result of such an addition could be the reason for the 
diversities in the observed circular polarization. Since our model was aimed at analyzing the 
intrinsic polarization properties of coherent radiation, we plan to consider the propagation 
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effects separately in tlie subsequent works. We speculate that the propagation origin of 
antisymmetric circular polarization is very unlikely but the symmetric circular polarization 
may be possible. 

6. Conclusion 

By taking into account of a detailed geometry of dipolar magnetic field lines, we have 
derived the polarization state of the coherent curvature radiation due to relativistic plasma 
in the pulsar magnetosphere, and drawn the following conclusions: 

1. We do confirm the previous results of Gil et al. (1990a,b; 1993) that coherent 
curvature radiation has basically antisymmetric type of circular polarization. Though 
the emission from a single bunch is highly polarized, the net emission from many 
bunches within the beaming region is less polarized due to the incoherent superposition 
of radiation fields. 

2. Based on the Stokes parameters of the curvature radiation we have deduced the 
polarization angle swing, i.e., the rotating vector model. 

3. Based on the coherent curvature radiation, we have achieved for the first time the 
result that the antisymmetric type of circular polarization is correlated with the 
polarization angle swing, and such correlations have been indeed found in the profiles 
of many pulsars. 

4. The addition of circular polarization with different signs and magnitudes at any given 
phase could be responsible for the wide diversity in circular polarization across the 
pulse. It is consistent with the earlier results (e.g., Gil et al. 1993, 1995). 

I thank J. L. Han and J. M. Rankin for illuminating discussions. 
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A.l. The Series Expansion Coefficients of Equation ( 1221) 



cosC {d,N,{9o) + d,N,{9o)) 
dj 

cosC {di {d, - /i) - 2d3d,N^ (60) - 2d2d,N, {Oq)) 

dl 

[di^df — Ad4^di2di — 2djdg,) cos( 

2df ' 
{dl (rfn - 3 (/4 + /s) - /e) - 2 (44 + Wis) dj + rf^) cosC 

c^iiVy (gp) - e,/3y (60) 
yo 72 



^x2 

^x3 

A. 

A. 

As 



yi 



zO 



A, 



z3 



^^At; (gp) - die^Py (60) + esATy (gp) + 6162 

rfi (eio - 2d,d^N; (gp)) + d, {2dl - d^d^) Ny (gp) + eg/^y (gp) 

rfi (3rfi (rfirfs - - 2d^) N'y (^p) + d^e^^ + 3ei4) - rfiCnATy (^p) + ei2/3y (^p 



sinC (-4A'x(go)-4A^z(^o)) 



A-2) 
A-3) 
A-4) 
A-5) 
A-6) 
A-7) 



A-8) 
A-9) 



. sinC (rfi(-rf3A^^(^p)-ci2A^a^o) + /i) + 24(4A'x(^o) + 4A^z(^o))) 

A^l — , (^A-iUJ 

_ (rfi(4rf4(rf9-/i) + 4(-c?io + /2 + 2/3) + 244)-6dl4)sinC 

ifr + dl (/e - rfii + 3 {f, + /s)) + 2c^f (44 + 34^0 - 34 (/2 + 2/3))) sinC 

(A-12) 



where 



dl = sin C /3x (6'o) + cos C /3z (^0) - 1 , 
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d2 = sinC -^x(^o) , (A-14) 

ds = /9.(^o)-cosC, (A-15) 

d, = sinC^;(^o)+cosC^^(^o) , (A-16) 

4 = sin C (^o) + cos C (^o) , (A-17) 

de = sinC/3x^'UM+cosC/3z^'UM , (A-18) 

d7 = di4 - 3dl , (A-19) 

= dsN^{eo) + d2N,{eo) , (A-20) 

= dsKi0o) + d2K{0o) , (A-21) 

rfio = dsN^ {Oo) + d2N^ (60) , (A-22) 

rfii = dsNj'^{0o) + d2N,^'Hdo) , (A-23) 

di2 = d9-iVz(^o)/3x(^o) + A^x(^o)/3z(^o) , (A-24) 

di3 = dio-27V^(^o)^x(^o) + 27V^(^o)^a^o)-A^z(^o)^x(^o)+ (A-25) 

di4 = ISdidid^ds - 2Adsdl - ddidjd^ , (A-26) 

ei = sinCiVx(^^o) + cosCiV,(^^o) , (A-27) 

62 = 2d4/3y (60) - di/3'^ {60) , (A-28) 

63 = di/3;(^o)-2rf4/3x(^o) , (A-29) 

64 = sinCA^^(^o) + cosCA^a^o) , (A-30) 

65 = cos C (^0) - 2^4 COS C (^o) + 2(^4 + 63 sin(C) , (A-31) 

66 = sinCA^x (^o) + cosCiV^'(^o) , (A-32) 

67 = sinCAx^^^ (^0) +cosCA"z^^U^o) , (A-33) 

68 = ciiiVy(3) (^0) - ei/3y(3) (^^) ^ (A_34) 

69 = 4^1^464 — 6(^461 + di (2^561 — diCe) , (A-35) 
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= di (diTV; (^o) - ei^;' (^o)) + (4^461 - 2die4) (^o) , 
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(A-46) 


f7 


= 24(^3^8 - 6di (^7 (/i - 4) + 3d4«) ■ 
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The expression of /3= {p^, py, P^} and the derivatives f3', f3", /3^^^ and (3^^\ which 
respectively represent the first, second, third, and fourth order differentiations with respect 
to 9 evaluated at 9o, are as follows: 

/3(^o) = {fi:(/iisin(r(^o)) + /i3cos(r(^o))), «:(/i2 sin(r(^o)) + /i4 cos(r(^o))), 

k(cos a cos(r(eo)) - h sin(r(^o))) } , (A-48) 

/3'(^o) = {^r'{eo){h,cosir{9o)) - h,smiri9o))), Kr'{9o)ih2cos{r{9o))- 

/i4sin(r(^o))), Kr'(^o)(-/i5cos(r(eo)) -cosa sin(r(^o)))} , (A-49) 

13" {9o) = \^n(T\9o)\h8sm(T{9o))-cos(P'{hgsm(T{9o))+sma cos(r(^o)))) + 
r"(^o)(/iicos(r(^o)) - /i6sin(r(^o)))), K{r"{9o){h2cos{r{9o)) - 
/i4sin(r(^o))) + r'(^o)'(/i2(-sin(r(^o)))-/i4cos(r(^o)))) , 
K{T'{9of{hsm{T{9o)) - cos a cos(r(^o))) - T"{9o){hcos{T{9o)) + 
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cosa sin(r(^o))))} , (A-50) 
/3(3)(^,) = |«;(r(3)(^o)(/ii cos(r(^o)) - /i6sin(r(^o))) + r'(^o)'(/i6sin(r(^o) - 

/i9cos(r(^o))cos0') + /i8cos(r(^o))) - {9o)r" {9o){hi sin(r(^o)) + 
/i3cos(r(^o)))), «:(r(3)(^o)(/i2cos(r(^o)) - /i4sin(r(^o))) + 

r'{0of{h^sm{r{0o)) - h2cos{r{0o))) - 3r'{0o)r"{0o){h2sm{r{0o)) + 
/i4cos(r(^o)))), /€(-r(=^)(^o)(/i5cos(r(^o)) + cos« sin(r(^o))) + 
r'(^o)'(/i5cos(r(^o)) + cosa sin(r(eo))) + 3r'(eo)r"(eo)(/i5sin(r(^o)) - 

COS a cos(r(^o))))} , (A-51) 

f3^'\0o) = [K{h, COS0' (sin(r(^o))(-3r"(^o)' + r'(^o)' - 4r(3)(^o)r'(^o)) + 
cos(r(^o))(r(^)(^o) - 6r'(^o)'r"(^o))) + /i8(sin(r(^o))(3r"(^o)' - 
r'(^o)' + 4r(3)(^o)r'(^o)) + cos(r(^o))(6r'(^o)'r"(^o) - r^^H^o))) + 
/i6(sm(r(^o))(6r'(eo)V'(eo)-r(^n^o)) + cos(r(^o))(-3r"(eo)' + 
r'(^o)' - 4r(3)(^o)r(^o)))), «:(/i2sin(r(^o))(-3r"(^o)' + r'(^o)' - 
4r(=^)(^o)r'(^o)) + /i2cos(r(^o))(r(^)(^o) - 6r'(^o)'r"(^o)) + 
/i4(sin(r(^^o))(6r'(^o)'r"(^^o) - r^^^^^o)) + cos(r(^o))(-3r"(^o)' + 
r'(eo)'-4r(3)(^o)r'(^o)))), «(cosa(sin(r(^o))(6r'(^o)'r"(^o) - 
r(^)(^o)) + cos(r(^o))(-3r"(^^o)' + r'(^^o)' - 4r^'\9o)r{9o))) - 
/i5(sin(r(^o))(-3r"(^o)' + r'(^o)' - 4r(3)(^o)r'(^o)) + 
cos(r(^o))(r(^)(^o) - 6r'(^o)'r"(^o))))} , (A-52) 

where 
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/i4 = sin a sin</)' , (A-56) 

/15 = sin a COS0 , (A-57) 

= sin a cos0' , (A-58) 

= sin a sin0' , (A-59) 

hs = sin0 sin0' , (A-60) 



(r/a \\ l + 3cos(2^o) 

cos(r(6'o)) = , , ^ , , (A-62) 



• f-pfa \\ 3 sin (2^0) . . ^ox 

sin(r(6'o)) = , (A-63) 



r (^o) = tan-i ( ^ ^^y. X ) ' (A-64) 



/ig = COS a COS0 , (A-61) 
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Note that the 0', which appears in the above equations, is just a variable for the rotation 
phase, and the prime (') on it does not represent any differentiation. 

By having known the derivatives of f3 from equations (]A-48p - (]A-52p . we can define the 
derivatives A?" evaluated at 6q : 
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A. 2. The series expansion coefficients of equation (I24h 



Co = gi (2(73 sin^ (^o) - 3^75 cos (^o) + 2(?2 + ^3 (log(2) - 2 log ((74))) , (A-73) 

ci = 3^1 ( sin (^0) (^3 sin (2^o) + 2^10 + 6 cos (20o)) + 

/€ cos (Fo) (sin (^o) - 3 sin (3^o)) ) , (A- 74) 

3c/i (^6 (1 + 3 cos (2^0)) + 2k cos (Fo) (5 + 3 cos (2^o)) (cos (^o) - 9 cos (3^o))) 
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4 (5 + 3 COS (2^0)) 

(A-75) 

if 4v/2 (28 sin (gp) + 9 (5 sin (3go) + sin {he,))) \ 
= 4^^ 1^-^^ + (5 + 3 cos (2^0))^/^ j ' ^^"''^ 



wliere 



cos (Fq) = sin a sm( coscp' +cosa cos^ , (A- 77) 

g2 = 6 + y31og(2 + Vs) , (A- 79) 

5^3 = 6k (sin C (sin sin 0' — COS a COS {/) COS 0' ) + sin a COS C COS (/) ) , (A-80) 

^4 = cos (^0) + ^5 + 3 cos (2^0) , (A-81) 

g5 = 4:K cos (Fo) sin^ (^0) + v^lO + 6 cos (2^0) , (A-82) 

ge = 7g3 sin (^o) + 3^3 sin (3^o) + 8^10 + 6 cos (2^o) cos (^o) , (A-83) 

= 9^3 cos (3^0) - 2/t cos (Fo) (sin (^0) - 27 sin (3^o)) • (A-84) 
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B.l. To Find Solution to Integrals in Equation (1251) 



Consider the integral 



By changing the variable of integration /i = (x/l) + m, and defining the constants / 
and m = —02/(803), we obtain 
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where Ai(z) is an entire Airy function of z with no branch cut discontinuities, and 
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where Ai'(2;) is the derivative of the Airy function Ai(2;). Therefore, we have 
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By differentiating equation (lB-21) on both sides with respect to ci we obtain 
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Differentiation of equation f lB-2|) on both sides with respect to C2 gives 
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Next, by differentiating equation ( 1B-2P on both sides with respect to C3 we obtain 
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Fig. 1. — Geometry for the calculation of radiation field at P, which is at a distance R from 
the source S. The magnetic axis rh is inclined with respect to rotation axis Qhj a. The sight 
line n impact angle with respect to rh is a. The colored curves represent the dipolar magnetic 
field lines plotted with r^, = 100 and azimuthal (0) increment of 30° for each field line, chosen 
rotation phase 0' = 0. The source position vector is r, velocity is v and acceleration is a. 
NS is the neutron star and C is an arbitrary field line. 
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Fig. 2. — Time t is plotted as a function of magnetic colatitude 9 of the bunch for different 
values of field line constant re- The normalization parameter P is the pulsar period. Given 
K = 1. 




Fig. 3. — Beaming regions specifying the range of magnetic colatitude 9 and azimuth at 
the three selected phases cf)' — —30°, 0° and 30°. The center of each region gives the values 
of 00 and ^o- Given a = 10°, /3 = 5° and 7 = 400. 
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Fig. 4. — Beaming regions specifying the range of magnetic colatitude 6 and azimuth 0. 
They are plotted for the full range of phase: -180° <(/)'< 180° with a step of 5°. The 
center of each region gives the values of 0o and 6*0- Given a = 10°, /3 = 5° and 7 = 400. 
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Fig. 5. — Simulated pulse profiles: in panels (a) and (b) intensity (/g), linear polarization (Lg) 
and circular polarization (K), and in lower panels (c) and (d) the corresponding polarization 
angle {ips) curves are plotted. Given P = 1 s and 7 = 400. Note that profiles are normalized 
with the peak intensity. 
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Fig. 6. — Simulated pulse profiles. Given P = 1 s and 7 = 400. = 0.1, 0p = 0°, and 
/o = 1 are used for the modulating Gaussian. 
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Fig. 7. — Simulated pulse profiles. Given P = 1 s and 7 = 400. For panels (a) and (c) 
= 1, 0p = 0°, and /o = 1, respectively, are used for the Gaussian. Similarly, for panels 
(b) and (d) a^f, = 1, (pp = 180, and fo = l are used, respectively. 
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Fig. 8. — Simulated pulse profiles. Given P = 1 s and 7 = 400. For panels (a) and (c) 
(j0 = 0.45, 0.32, (pp = 0°, ± 60°, and /o = 1, 0.9, respectively, are used for the Gaussians. 
Similarly, for panels (b) and (d) = 0.45, 0.32, 0p = 180°, 180° ± 60°, and /o = 1, 0.9, 
respectively, are used from the central component to the outer one. 
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Fig. 9. — Simulated pulse profiles. Given P = 1 s and 7 = 400. For panels (a) and (c) 
= 0.14, 0.10, 0.07, 0p = 0°, ± 30°, ± 50°, and /o = 1, 0.75, 0.5, respectively, are 
used for the Gaussians. Similarly, for panels (b) and (d) cr^ = 0.14, 0.06, 0.03, (pp = 
180°, 180° ± 16°, 180° ± 26°, and /o = 1, 0.75, 0.5, respectively, are used from the central 
component to the outermost one. 
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Fig. 10. — Simulated pulse profiles. Given P = 1 s and 7 = 400. For panels (a) and (c) 
= 1, 0.1, 0.05, (pp = 0°, ± 65°, ± 75°, and /o = 1, 0.9, 0.8, respectively, are used 
for the Gaussians. Similarly, for panels (b) and (d) = 1, 0.1, 0.05, (pp = 180°, 180° ± 
65°, 180° ± 75°, and /o = 1, 0.9, 0.8, respectively, are used from the central component to 
the outermost one. 



